Role of pulling direction in understanding the energy landscape of proteins 
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Single molecule force spectroscopy provide details of the underlying energy surfaces of proteins 
which are essential to the understanding of their unfolding process. Recently, it has been observed 
experimentally that by pulling proteins in different directions relative to their secondary structure, 
one can gain a better understanding of the shape of the energy landscape. We consider simple lattice 
models which are anisotropic in nature to study the response of a force in unfolding of a polymer. 
Our analytical solution of the model, supported by extensive numerical calculations, reveal that 
the force temperature diagrams are very different depending on the direction of the applied force. 
We find that either unzipping or shearing kind transitions dominate the dynamics of the unfolding 
process depending solely on the direction of the applied force. 



PACS numbers: 64.90.+b,36.20.Ey,82.35. Jk,87.14.Gg 



I. INTRODUCTION 

Recent experiments have shown that the energy land- 
scapes governing the mechanical unfolding of proteins are 
different from the energy landscapes governing the ther- 
mal or chemical unfolding P, H, S] • Until recently it 
was believed that all residues of any given protein con- 
tribute equally to the thermodynamic stability of glob- 
ular proteins, with the stability of the protein structure 
being provided by its hydrophobic core @, 0). It was 
therefore expected that the mechanical unfolding of a 
protein should be insensitive to the direction of the ap- 
plied force. Surprisingly, recent experiments have shown 
that a protein's resistance to unfolding depends strongly 
on the pulling direction || 0, [l(| [H[ ■ Varying the pulling 
direction also gives Angstrom precise information about 
the structure of biomolecules [H, EH ■ 

Force extension curves obtained by single molecule ex- 
periments (unfolding of proteins, stretching of DNA) are 
often described by the worm like chain model (WLC 
model) or freely jointed chain model (FJC model) [3, 
[l5j . The mechanical properties of the FJC and WLC 
models are well understood, and the usefulness of these 
models stems from their simplicity, which means that the 
analytic expression for the extension can be written in a 
very simple form [lril [l7j . The only parameter appear- 
ing in the description of these models is the persistence 
length of the polymer chain. The force-extension curves 
obtained from these models [l7l [IH show an excel- 
lent agreement with the main features observed in early 
experiments [J H i, H . Since the FJC and WLC models 
are isotropic, in order to study the resistance of a pro- 



tein to unfolding when the pulling geometry is changed, 
one would have to use different persistence lengths (along 
different directions) to fit the force extension curves 
[l,S[i3- S° when a polymer has been pulled along the x- 
direction one may use a certain persistence length which 
reproduces the observed force-extension curve. If the 
pulling direction is changed another persistence length is 
used to fit the new force-extension curve. However, the 
underlying model is still isotropic and the procedure is 
somewhat ad hoc and does not provide a unified approach 
to the modeling of the anisotropic behavior [ijj [2(| seen 
in experiments. Equally importantly these simple models 
do not incorporate excluded volume effects in its descrip- 
tion. 

In a recent paper Kumar and Giri [2l[ used par- 
tially directed self-avoiding walks (PDSAWs) to model 
anisotropic biomolecules. When a force is applied to one 
end of the chain it undergoes a transition from a folded 
to an unfolded state. Numerical studies based on a chain 
size of N — 30 showed that when the chain is being 
pulled along the preferred direction of the polymer, the 
nature of the transition is akin to unzipping, but when 
force is applied perpendicular to the preferred direction, 
the transition is akin to shearing. 

In section II we provide the analytical solution of the 
model and show that a change in the pulling direction 
gives rise to different phase diagrams even in the ther- 
modynamic (infinite chain length) limit. In section III 
we report on an algorithmic breakthrough which has en- 
abled us to increase the chain length to N = 150 which 
is five times longer than in the previous study [21j . The 
longer series are used to obtain numerical estimates for 
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the phase boundary in excellent agreement with the ex- 
act results. In this section we also discuss the relevance 
of two ensembles namely, the constant force ensemble 
(CFE) and the constant distance ensemble (CDE) which 
are appropriate when describing different experimental 
setups. With the enhanced information about the exact 
density of states of longer chains, we are able to study 
precisely finite-size effects which are crucial to all sin- 
gle molecule experiments. In particular we find marked 
differences between the cases where the force is applied 
along the x- or y-directions, respectively, and we explain 
these observed differences using simple heuristic argu- 
ments. The PDSAW is a very restricted model and in 
order to further examine the role of anisotropy we in- 
troduce an anisotropic self-attracting self-avoiding walk 
(ASASAW) model and using exhaustive exact enumer- 
ation data for chains up to length N = 48 study the 
combined effect of anisotropy and pulling in section IV. 
The paper ends with a brief discussion in section V. 



II. PDSAW MODEL AND PHASE DIAGRAM 
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FIG. 1: (Color online) Schematic representation of a PDSAW 
on the square lattice with one end kept fixed and the other 
end subjected to a pulling force along (f x ) and perpendicu- 
lar (f y ) to the preferred direction. At low temperature the 
conformation mimics a /3-sheet. 



The PDSAW is a self-avoiding walk in which steps 
along some direction (say, the negative x-direction) are 
forbidden. We study this model on a two dimensional 
square lattice. An attractive energy e u is associated 
with each non-bonded nearest neighbor [22]. At very 
low temperatures the attractive interactions dominate, 
and the PDSAW is in a collapsed phase, where the den- 
sity of monomers in the bulk is close to one. The typical 
configurations of the chain mimics the /3-sheet [19j (see 
Fig. [I}. At high temperatures thermal fluctuations are 
strong enough to break some of the interaction bonds and 
with increasing entropic contributions to the free energy 
the chain can be in an extended phase. A force F may 
be applied along some fixed direction giving rise to a 
stretching energy E s arising from the applied force F 



E, 



-F ■ a 



(1) 



where a is the end-to-end vector. In this paper we con- 
sider only the cases where the force is applied either along 
the x-direction (f x ) or the y-direction (f y ). 

The phase boundary between the collapsed and the ex- 
tended phase can be obtained by calculating the macro- 
scopic shape of the collapsed phase at low tempera- 
tures. Let uj x — exp(/a;/r), uj y = exp(/j,/T) and 
u = exp(e u /T), where T is the temperature. We set the 
Boltzmann constant ks equal to I. The energy of the 
configuration which mimics the /3-sheet is (see Fig. [T]) 

b-2 

E = -e u N + ^( 0l + a 2 + 2b) + ^ £ \y j+2 - yj\, (2) 

j=o 

where, as illustrated in Fig.Q] cii = \yi — yo| is the height 
(number of steps in the y-direction) of the first column, 



a-2 = \Vb -2/6-1 1 is the heig 
the number of steps in the 
Following the approach 
the sum appearing in Eq. 
The macroscopic shape of 
the value of lu x at which 
unstable gives the critical value, u x , in terms of u> y and 



ht of the last column and b is 
x-direction. 

developed in [23|, we express 
@ as a surface tension term, 
the polymer is calculated and 
the collapsed phase becomes 



(u-1) 



[(Vu - ^){y/u y JUJ v ' - I)] . (3) 



Eq. ([3]) gives the complete three dimensional phase 
boundary. For f y = Q (uj y = 1) Eq. (JU) reduces to the 
following simple expression for the critical force when 
pulling along the a;-direction 



(4) 



This coincides with the result obtained by transfer ma- 
trix methods [2(|. However, that method is not easily 
generalized to the case where force is applied along the 
y-direction. Eq. ([3]), on the other hand, is sufficiently 
general to yield an expression for the critical force when 
pulling along the y-direction. Setting f x = (or lj x = 1) 
yields the critical force when pulling along the y-axis: 



1 + u 2 + \/T+ 2u 2 - 4u 3 + u 4 
2w 3 / 2 



(5) 



For zero force (i.e. ui x = u> y = I) the expression re- 
duces to the exact value u c = 3.3829 . . . (T = 0.8205 . . . 
if we set e u = 1 ) at which thermal unfolding occurs. The 
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FIG. 2: (Color online) The globule-coil phase boundary in the 
force-temperature plane. The response to the force when ap- 
plied along the y-direction (squares) is distinctly different to 
the case where the force is applied along the i-direction (cir- 
cles). The continuous lines show the exact phase boundaries 
which are in excellent agreement with our numerical results 
for finite N = 150. 



exact phase boundaries obtained from Eqs. Q and ([5|) 
are shown in Fig. [2] and compared to numerical results 
obtained for N = 150. 



III. NUMERICAL ANALYSIS 
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FIG. 3: Top: Fluctuations in the number of nearest neigh- 
bor contacts (x) as a function of temperature for different N. 
Inset shows how the peak height increases with N at the tran- 
sition temperature; Bottom: The variation of the transition 
temperature as TV — » oo. The dot on the y-axis marks the 
exact transition temperature. 



Our numerical studies of the process of force induced 
protein unfolding are based on the exact enumeration of 
all possible conformations of the chain. This approach 
not only enables us to obtain precise estimates for the 
phase boundary, but also allows us to study the interme- 
diate states crucial to events occurr ing during the me- 
chanical stretching process [HI, US Ha]- In this study 
we have combined the power of parallel processing and 
transfer matrix calculations to extend the limit of earlier 
studies not just by one or two monomers but five times 
over! Thus we have extended from the previous longest 
chain length 2l| N = 30 to much longer chains of up to 
N = 150 steps. The partition function is 

Z N = C N {N p , ax , ay )u N *^uj«y, (6) 

Np,a x ,cty 

where Cn is the total number of /V-stepped PDSAWs 
having N p non-bonded nearest neighbor pairs with end- 
to-end vector (a x ,Oy). 

The value of the transition temperature (for fixed force 
and finite N) can be obtained from the fluctuations in 
the number of non-bonded nearest neighbor interactions. 
The fluctuations are defined as \ — (Np) — (N p ) 2 , with 



the fc'th moment given by 



E 



C N (N pi a x , a v )u N PUx x cOy v 



(7) 



When the fluctuations \ are plotted as a function of tem- 
perature the resulting graph (Fig. [3^) has a peak at the 
transition temperature. The inset shows that the height 
of the peak in the fluctuation curve grows as a power-law 
with N, this being the hall mark of a phase transition. 
By setting e u — 1 we obtain (for N = 150) at zero force 
(u) x = u) y = 1) the transition temperature T = 0.8446. 
This is already quite close to the value T = 0.8205 ob- 
tained in the thermodynamic limit (TV — ► oo). The bot- 
tom panel in Fig. [3] shows the variation in the transi- 
tion temperature plotted against 1 /N and we see a pro- 
nounced curvature for large N making an extrapolation 
to the limit N — > oo very difficult. 

In Fig. [5] we show the force-temperature diagrams for 
N = 150 when pulling along the x- and y-axis (circles 
and squares, respectively). These numerical estimates 
for the phase-boundary are in excellent agreement with 
the analytical results (shown as solid lines) obtained from 
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FIG. 4: Plots of the average extension (x)/N {{y)/N) vs the applied force F for chain lengths N = 150, 100 and 50 at fixed 
temperature T = 0.01 with the force along the a;-direction (y-direction) . 



Eqs. ((3]) and ([5]). Thus exact enumeration technique 
works quite well even for small N as far as the phase 

diagram is concerned. This has also been seen earlier (a) I I I I I I I I I I I I (d) 

[13] in case of the surface adsorption. 



A. Constant Force and Constant Distance 
Ensembles 

Many biological reactions involve large conformational 
changes providing a well denned mechanical reaction co- 
ordinate (e.g. the end-to-end distance of a bio-polymer), 
which has been used in force spectroscopy experiments 
to follow the progress of the reaction [l], lf|. Within the 
constant force ensemble (CFE) the theoretical predic- 
tions based on exact enumeration results provide qual- 
itative description of the outcomes of such experiments 
[2ll [28[ (like the existence of intermediate states and dif- 
ferent pathways to unfolding) . In Fig. 0] we have plotted 
the average extension as a function of the applied force. 
Fig. fjji shows the result when the force is applied in the 
x-direction and Fig. [5Jd when the force is applied in the 
y-direction. In both cases we see the existence of several 
plateaus showing that there are many intermediate states 
between the collapsed and fully extended state. A major 
difference is that the number of these intermediate states 
seems to grow much faster with the chain length N when 
pulling in the x-direction as compared to pulling in the 
y-direction. 

The appearance of the plateaus in the force-extension 
curves as well as the differences between the two cases 
have a simple heuristic explanation. At very low temper- 
ature there are no entropic contributions to the free en- 
ergy which is thus dominated by the energetic terms aris- 
ing from nearest-neighbor interactions and the applied 
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FIG. 5: Schematic diagrams of the configurations giving rise 
to the plateaus in the force-extension curves when pulling in 
the x-direction (figures (a)-(c)) and along y-direction (figures 

(d)-(f))- 



force. When pulling in the x-direction we have to maxi- 
mize the quantity G = Fx + N p , where N p is the number 
of contacts. For F — ► this means just maximizing m, 
and this is achieved by zigzag conformations inscribed 
in a rectangle which is as close to a square as possible 
so that (a;) /N is close to \/N. As the force is increased 
other configurations may become energetically favored, 
essentially corresponding to decreasing the width of the 
rectangle by a unit and increasing the length correspond- 
ingly (a more detailed explanation can be found in [29| ; 
see also (30[). So the plateaus seen in Fig. [4^, arise from 
inscribing a PDSAW of fixed length N in rectangles, 
going from right to left, of width w — 2,3, .. . y/N as 
evidenced by the fact that these plateaus have height 
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FIG. 6: Plot of the average force (F) vs the elongation at different temperatures T for chain length N = 150. The marked 
difference due to a change in pulling direction is clear from these plots. 



(x)/N = 1/w = 1/2,1/3,.... The types of confor- 
mations giving rise to the right-most plateaus are rep- 
resented schematically in Figs. [5^,-c. When pulling in 
the y-direction we again get a competition between the 
force which tends to maximize the distance of the end- 
point in the y-direction and interactions which try to 
inscribe the PDSAW into a square. As illustrated in 
Figs. [5ji-f the conformations giving rise to the right- 
most plateaus are now those where the PDSAW is in- 
scribed in a rectangle of length I = 3, 5, 7, . . . , y/~N. This 
can be seen in Fig. |4p where these plateaus have height 
{y)/N = l/l = 1/3,1/5,.... Even lengths cannot oc- 
cur because for these the end-point would have a y- 
component close to giving only a small energetic con- 
tribution from the applied force. The upshot is that for a 
given length the number of plateaus when pulling in the 
x-direction will be about twice the number of plateaus 
when pulling in the y-direction. 

Our longer series data also allows us to go a step fur- 
ther in analyzing the model in the constant distance en- 
semble (CDE). This is the ensemble best suited to the 
analysis of experiments performed using apparatus such 
as atomic force microscopes. In Fig.|B]we have plotted the 
average force (F) [28| as a function of the extension when 
pulling in either the x- or y-direction. Striking differences 
between the two cases are obvious from these plots. Most 
obviously we note that at the low temperature T = 0.1 
the force extension curve obtained for pulling in the x- 
direction shows unzipping like transition characterized 
by having smooth plateaus. However, when the chain 
is being pulled in the y-direction (at the same tempera- 
ture) , the force extension curve exhibits "saw-tooth" like 
oscillations indicating that the transition is akin to shear- 
ing. This clearly demonstrates that even the nature of 
the unfolding transition can change dramatically solely 



by varying the pulling direction. 



IV. ANISOTROPIC-SELF-ATTRACTING-SELF- 
AVOIDING 
WALKS 

The model described above has an exact solution and 
gives qualitative features similar to experiments. How- 
ever, the physical constraints imposed by experimental 
set-ups has not been taken fully into account. For ex- 
ample, in experiments using atomic force microscopes, 
receptor and ligand molecules are attached to a sub- 
strate and a transducer, respectively, which introduce 
anisotropy in the systems. Also steric constraints due to 
the confinement imposed by the experimental set-up can 
lead to a loss of entropy and thus result in modifications 
to the behavior of the chain. Furthermore the constraint 
of partially directing the walk so it cannot take steps 
in the negative x-direction is not really appropriate to 
most experimental conditions. Such a constraint is typi- 
cally only appropriate when the system is under flow or 
constant external field [U [HJ . A more realistic model 
of polymers is provided by self-attracting-self avoiding 
walks (SASAWs), which is well suited to the modeling of 
a linear chain in a poor solvent [HI, However, this 
model is isotropic in nature. In this paper we introduce 
anisotropy to the model. We achieve this by consider- 
ing nearest neighbor interactions with different strengths 
along the x- and y-directions as shown in Fig. [7^. This 
is in accordance with real proteins where the interactions 
along different directions can differ by orders of magni- 
tude [III] . Hence by changing the strength of the interac- 
tions one can vary the degree of anisotropy in the system 
and study the effect this has on the force-temperature 
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FIG. 7: (Color online) Schematic representation of ASASAWs 
on the square lattice, (a) Different nearest neighbor interac- 
tions along the x- and y- direct ions introduce anisotropy in 
the model polymer, (b) represents ground state conforma- 
tions which dominate the partition function similar to Fig. [T] 
when subjected to force along the ^-direction, (c) Same as 
(b) but resembles Fig. Q] when subjected to force along the 
y- direct ion. 



phase diagram. The partition function for anisotropic- 
self-attracting-self-avoiding walks (ASASAWs) may be 
written as 



N 



(N p , 



E 



C N (N 1 



px i N py i 



\a\)% 



(8) 

Here Cn is the total number of ASASAWs of N steps 
having N px and N py nearest neighbor pairs along the 
x- and y-directions, respectively, while u x = exp[— /3e x ] 
(u y = exp[— 0e y ]) are the Boltzmann weights associ- 
ated with the nearest neighbor interactions between non- 
bonded monomers along, respectively, the x- and y- 
directions. 

The qualitative as well as quantitative features of the 
force-temperature phase diagram as shown in Fig. [8] for 
the isotropic case (i.e . u x — u y — u) remain the same as 
those reported in [28| • To make our model closely resem- 
ble the PDSAW discussed above, we set either u x — 1 
or u y = 1 and apply a force along the x-direction. For 
u y = 1 the ground state is dominated by configurations 
similar to the one shown in Fig. [7](b) . This closely re- 
sembles the PDSAW model (Fig. Q} subjected to a force 
along the x-axis and the nature of the transition is akin 
to an unzipping transition. On the other hand for u x = 1 



N=48 (ASASAW) 
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FIG. 8: (Color online) The force-temperature phase diagram 
for ASASAWs. The solid line represents the phase boundary 
for the isotropic SASAWs. 



the ground state is dominated by configurations like the 
one shown in Fig. [7^c). This closely resembles the con- 
figurations of Fig. [T] with a force along the y-axis and 
consequently the nature of the transition is akin to a 
shearing type transition. It should be noted that the 
number of ASASAWs, Cat, is different from the number 
of PDSAWs even if we put u y = 1. We observe a sig- 
nificant decrease in the transition temperature compared 
to SASAWs. This is caused by a decrease in entropy 
which allows the polymer to acquire conformations sim- 
ilar to a /3-sheet. The effect of surface confinement is 
also evident from these plots. The transition temper- 
ature is found to be less for u x = 1 in comparison to 
u y = 1. Another major difference between the ASASAW 
and SASAW phase-diagrams is that there appears to be 
no re-entrance in the ASASAW model when u x = 1 or 
u y = 1. It is well known that the re-entrant behavior of 
the SASAW model is due to a non-zero ground state en- 
tropy, that is, the number of configurations in the ground 
state grow exponentially with N (these configurations are 
just Hamiltonian walks inscribed in a rectangle with side 
lengths as close as possible to V~N). For the ASASAWs 
we study here the ground-states are, as argued above, 
folded sheets as in Fig|7fb) and (c), respectively. The 
number of such configurations do not grow exponentially 
with TV and are thus too few to give a non-zero ground 
state entropy, which is why we do not see re-entrance. 
However, if both u x and u v are greater than 1 then one 
may see a re-entrant phase diagram. 

The force extension curve shown in Fig. [9] not only 
shows the existence of intermediate states but also the 
emergence of new states depending on the direction of the 
interactions in accordance with the findings for PDSAWs. 
As the temperature is increased we observe in Fig. 10 (as 
expected) that the intermediate states are washed out 
due to the resulting increase in entropy. We also observe 
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FIG. 9: Left panel: The average scaled elongation {x)/N as a function of the force with interactions along the a;-direction at 
fixed low temperature T = 0.01 for different chain lengths; Right panel: Same as the left panel but for interactions along the 
y- direct ion. 
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FIG. 10: Left panel: The average scaled extension (x) /N 
along the x-direction as a function of the force at different 
temperatures with chain length N = 48; Right panel: Same 
as the left panel but with interactions along the y-direction. 




FIG. 11: Left panel: The average force (F) vs the extension 
x with interactions along the a;-direction alone at different 
temperatures with chain length TV = 48; Right panel: Same 
as the left panel but with interactions along the y-direction. 



that the system attains a higher stability against the force 
when the interactions are along the y-direction and that 
intermediate states survive even at high temperatures as 
may be seen from Fig. [TU1 

Finally we analyze the ASASAW model in the con- 
stant distance ensemble. In Fig. [TT] we have plotted 
the average force (F) [28| as a function of the extension 
when pulling in the x-direction with interactions either 
along the x- or y-direction. The results are qualitatively 
very similar to those obtained for the PDSAW model 
(see Fig. In particular we note that the oscillation 
at low-temperature are much more pronounced in the 
case where the interactions are in the y-direction. As ex- 
plained above (see Fig. [7]) there is a close correspondence 
between the ground-state (low-temperature) configura- 
tions of the ASASAW with interactions in the x- or y- 
directions and the PDSAW with a force applied along the 
x- or y-direction, respectively. This close correspondence 
is reflected in the similarity of the results. 



V. CONCLUSION 

In summary we have demonstrated via the analytical 
solution for the phase boundary of PDSAWs as well as 
through high precision numerical calculations that chang- 
ing the pulling direction can change the nature of the un- 
folding transitions and give rise to very different phase di- 
agrams. The emergence of many new intermediate states 
suggests that there may be many pathways to the unfold- 
ing of a polymer. The findings for the PDSAW model 
have been confirmed by our study of ASASAWs. The 
inclusion of anisotropy in the traditional SASAW model 
of polymers gives unequivocal evidence that the features 
observed here for the PDSAW model are in fact generally 
true and not an artifact of the PDSAW model. In future 
additional work like extensive Monte Carlo simulations 
and molecular dynamics studies are needed to better un- 
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derstand the role played by anisotropy and how varying 
the pulling direction can help us gain a better and deeper 
understanding of biomolecules. 
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